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NONLINEAR CONTROL SYSTEM DESIGN METHODS 


Consider some of the control system design methods for plants with nonlinear 
dynamics. If the nonlinearity is weak relative to the size of the operating region, 
then the linear methods apply directly. Fixed-gain design may be feasible even for 
significant nonlinearities. It may be possible to find a single gain which provides 
adequate control of the linear models at several perturbation points. If the non- 
linearity is restricted to a sector, that fact may be used to obtain a fixed-gain 
controller. Otherwise, a gain may have to be associated with each perturbation point 
p^. A gain schedule K(p(v)) is obtained by connecting the perturbation points by a 
function, say p(v), of the scheduling parameter v (i.e., speed). When the sche- 
duling parameter must be multidimensional, this approach is difficult; the objective 
of our research is to develop an easier procedure. 


X = f (x, u, z) 

1. SMALL OPERATING REGION - LINEAR METHODS 

P = Kr u o } ' S * = (fj) 5x + (f^) 5u, 5u = K(p)5x 
P P 

2. FIXED-GAIN DESIGN 

a. ONE GAIN FOR ALL {pj 

b. ONE GAIN FOR SECTOR NONLINEARITIES 

3. GAIN SCHEDULING 

Pj = p(vj), K(p(v)), v = SCHEDULING PARAMETER 

4. NONLINEAR TRANSFORMATION OF STATE AND 
CONTROL 
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SIMPLIFICATION THROUGH COORDINATE CHANGE 


We attempt to simplify the design problem by simplifying the representation of 
the plant. Here is an example. The state equation is nonsingular, but a fixed-gain 
design is impossible. The nonsingular matrix E represents rotation in the plane 
of x through the angle <J>. The change of control coordinates from u to v results 
in a globally constant, linear system. Even an approximate cancellation of E could 
be quite helpful, since the resulting system may be nearly constant and the fixed- 
gain design may be applicable. 


xeR^, ueR^ EE r = I 
x = E (<//{t)) u 

$ = O — ► x = u, k = K q x, STABLE 
4/ = 7T — ► x = -u, x= -K q x, UNSTABLE 

NEW CONTROL COORDINATES: 

v = E{\p(t)) u ► k= v, GLOBALLY 

APPROXIMATE TRANSFORMATION: 

v = E (\p) u ► x = E (5i//(t)) v FIXED GAIN MAY WORK 
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BRUNOVSKY CANONICAL FORM 


The first step in the design approach is to try to transform the given system into 
something more simple, i.e., ideally, a set of decoupled strings of integrators. This 
set is called the Brunovsky canonical form for controllable, constant, linear systems. 
The n um ber of strings is given by the number of control axes m, and the number of 
integrators (dots in the figure) is given by a Kronecker index K^. In general, 

= n. Two examples are shown for the case of n = 12 and m = 4. According to 
linear theory, the set of Kronecker indexes is invariant under nonsingular transforma- 
tions and feedback. 


GIVEN x = Ax + Bu, xeR n , ueR m , (A, B) CONTROLLABLE 

CAN BE TRANSFORMED WITH 

y = Tx , v = Rx + Wu , (T, W) NONSINGULAR 

INTO m DECOUPLED STRINGS OF INTEGRATORS , LENGTH = KRONECKER INDEX 


EXAMPLE: n=12, m=4 


v-j *»• • • • 


v-j — fr 


v 2 


** • • — + 


v 2 — • • 


v 3 • 


v 3 — • • 


V4 • 


V4 — *-♦ # • 


Kl = (4,4,2, 2) 


Kl = (3,3,3,3) 
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OUTLINE OF DESIGN PROCEDURE 


In principle, the design procedure is to transform the natural representation 
of the plant into the corresponding Brunovsky form, then design a control law 
v = g(y) for the canonical system, and finally pull the law back into the natural 
coordinates in terms of which law must be implemented. 


NATURAL COORDINATES 
x = f(x, u) 


u = W(x, v) 
y = T(x) 


CANONIC COORDINATES 
y = A 0 y + B 0 v 


v = y 5 


« • • • 


• • • 


« • • • 


• 

Kl = (4,4,4,2) 


u = W(x.g(T(x))) 


v = g(y) 
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DEVELOPMENT OF THE DESIGN APPROACH 


Before this design approach becomes practical, several issues must be resolved. 

(1) When can a given system be transformed into a linear model, 
how does one construct the required transformation, and how 
feasible is it to implement the resulting algorithm on flight 
computers? 

(2) What is the structure of the complete control system that 
includes the linearization step? 

(3) How robust is the resulting design? 

(4) How can the constraints on control and state be enforced? 

These issues have been explored both theoretically and experimentally. The results 
are summarized in the following figures. 


1. TRANSFORMATION 

a. EXISTENCE 

b. COMPUTATION 

c. IMPLEMENTATION 

2. STRUCTURE OF COMPLETE CONTROL SYSTEM 

3. COMPLEXITY AND ACCURACY OF MODEL - ROBUSTNESS 

a. EXACT STATE SPACE 

b. TRUNCATED STATE SPACE 

4. ENFORCEMENT OF DESIGN CONSTRAINTS 

5. EXPERIMENTS 

a. FORTRAN 

b. REAL-TIME ( FLIGHT COMPUTER, HYDRAULICS) 
SIMULATION 

c. FLIGHT 
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LIE BRACKETS 


The key theoretical result from the existence of linearizing transformations is 
the result of work by Krener (ref. 1), Brockett (ref. 2), Jakubcyzk and Respondek 
(ref. 3), and Hunt and et al. (ref. 4). The necessary and sufficient conditions are 
best expressed in terms of lie brackets. A lie bracket (f , g) constructs a new vector 
field from the old ones f and g. A set is involutive if the brackets do not create 
new directions. 


VECTOR FIELDS 

f, g: R n -*■ R n 



LIE BRACKETS 

[f, gl = |&f- 

0X 

3f 

^-g 

9x 


INVOLUTIVE SET 

l h 1 h rl 

IF 



[hj.hj] G SPAN | h-j, . . 

...h r } 
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CONDITIONS FOR LINEARIZABILITY 


1. It is necessary to find control coordinates which appear linearly in the state 
equation. For aircraft this means angular acceleration instead of ailerons, 
elevator, and rudder. 

2. The resulting system must have linear-like controllability. 

3. The fields {f, g]_» § 2 " * ** must satisfy an involutivity condition. For 

example, let n = 6, m = 2, and the Kronecker index set KI = (3, 3). Then, the 
six vector fields must span the state space and the first four must be involutive 


1. x = f(x, u) — ► x = f (x) + 2 g t (x) u* 

2. LINEAR - CONTROLLABLE 

3. INVOLUTIVE 


EXAMPLE: x = f{x)+g 1 (x) u.,* + g 2 (x) u 2 * 

KI = (3, 3) 



6 DIMENSIONAL 
SPAN 
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STRUCTURE OF THE COMPLETE CONTROL SYSTEM 


The structure has the form of an exact model follower. The linearizing transfor- 
mation is done in the WT-map. Through this map the plant is seen as a set of 
decoupled strings of integrators. The same set is employed as the dynamics (A Q , B 0 ) 
of the model servo, where the desired motion is defined by means of the input r* and 
the, in general, nonlinear control law. The control of modeling inaccuracies and other 
disturbances is accomplished by the regulator which operates on the error ey and 
outputs corrective control 6v. It may be noted that the regulator works into the 
simple canonical dynamics and, in effect, the gain scheduling is done automatically 
by the WT-map. 



MODEL SERVO REGULATOR WT-MAP AIRCRAFT 
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HELICOPTER EXPERIMENT 


The objective of this experiment is to investigate the effectiveness and 
realism of the design approach and to uncover potential problem areas. The current 
work is with the UH1H helicopter equipped with the VSTOLAND avionics system including 
the Sperry 1819B flight computer. The model used in the design is a rigid-body 
nonlinear force (f^), and moment (f^) generation process. Inertial coordinates 
(r, v) of position and velocity vectors, body attitude matrix C, and angular veloc- 
ity 03 form the state. The moment controls u^ 1 
the pedals. The collective is the power control 


are the roll and pitch cyclic and 
u^ , 


STATE: 


x =| V 1 e X C R 3 X R 3 X SO(3) X R 3 
C 


VCO/ 


CONTROL: 


(u 3 

u=lp 1 e U C R 3 X R 


STATE EQUATION: 


r = v 

V = f F (x, u) 
C = S(co) C 

CO = f^(x, u) 
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A CANONICAL MODEL FOR THE HELICOPTER 


The canonical model chosen is shown in the figure. There is a pair of strings, 
each four integrators long. This pair represents the two horizontal channels 
(r x , r y ) . In addition, there is a two-integrator string for altitude h. The 
fourth string is for the heading 4 >. The canonical controls are the second deriva- 
tives of horizontal acceleration, vertical acceleration, and yaw acceleration. The 
transformation is computed by means of two Newton-Raphson. trim .routines . First, the 
controls are computed and yield the given accelerations (w^, h) at the given state. 
Then the attitude C is computed for the given horizontal acceleration. The 
Jacobian matrices needed by Newton-Raphson are computed numerically. 


CANONICAL MODEL 


h 

& 


a x a x v x r x 

** • • • 


& — • 4 * 

Kl = (4,4,2, 2) 


TRANSFORMATION 


a. ON-LINE NEWTON-RAPHSON MOMENT TRIM 

cM < -- « ' - ‘ ' 


f (r,v, C, w b , u) = w bc 

(r,v,C,w b ,u) = h c | 
b. ON-LINE N-R FORCE TRIM 


-►u c = h M (v, C, w b , h c , w bc ) 


( r, v, C, o, u_) = at. 


(0, d) = h F (r, v, t// , a h , o, o) 
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SUMMARY - UH1 EXPERIMENT 


The experiment has progressed to the point that the following observations can 
be made. 


(1) The implementation of the control scheme is practical. The 
complete code takes less than 22 msec on the Sperry 1819B 
flight computer 

(2) Both FORTRAN and real-time simulations (real flight computer 
and hydraulics) are consistent with theory. The tracking 
accuracy along the trajectory including hover, climb, descent 
and high-speed flight indicates that the UH1 model routinely 
used for manned simulations is linearizable to a degree 
where a fixed-gain controller is possible. There is robust- 
ness with respect to weight, center of mass, moment of iner- 
tia, and force and moment models, but there is sensitivity 

to errors in wind estimates. The actuator dynamics (15 
rad/sec) may be commuted with the TW-map for regulator 
bandwidth below 2 rad/ sec 

(3) Further research is needed to develop rigorous methods for 
including high-frequency dynamics and explicit enforcement 
of control and state constraints 


1. IMPLEMENTATION IS PRACTICAL 

SAMPLING TIME = 50 msec 

NEW JACOBIANS FOR N - R TRIM EVERY 5TH SAMPLE 
COMPLETE CODE = 22 msec ON SPERRY 1819B 

2. FORTRAN AND MANNED SIMULATIONS ARE CONSISTENT WITH 
THEORY 

TRACKING ACCURACY IMPLIES L IN EAR IZABILITY 

ROBUSTNESS 

ACTUATOR DYNAMICS 

3. PROBLEM AREAS 

METHOD FOR INCLUDING SERVO DYNAMICS 

METHOD FOR EXPLICIT ENFORCEMENT OF CONSTRAINTS 
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